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^^ . We analyze the anomalies of superconducting state (both s and d-wave pairing) in the model 

fy^ ' of the pseudogap state induced by Heisenberg spin fluctuations of antiferromagnetic short - range 

^3 I order, and based on the scenario of strong scattering near the "hot spots" on the Fermi surface. We 

Cn , present microscopic derivation of Ginzburg - Landau expansion, taking into account all Feynman 

r^ ' graphs of perturbation theory for electron interaction with fluctuations of short - range order and 

I— { ' in the "ladder" approximation for electron scattering by normal (nonmagnetic) impurities. We 

determine the dependence of superconducting critical temperature Tc and other characteristics of a 
superconductor on the parameters of the pseudogap and impurity scattering. It is shown that within 
this model it is possible to explain the typical phase diagram of high - temperature superconductors. 
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;Ll ■ INTRODUCTION 

^ ■ 

c/j ' One of the main problems in the physics of high ~ temperature superconducting cuprates remains the theoretical 

understanding of typical phase diagram of these compounds JJ. Especially interesting is the clarification of the 

2 ' nature of the pseudogap state which is observed in a wide region of temperatures and concentration of carriers J3] , 

and is obviously crucial for the formation of electronic properties both in normal and superconducting states. Despite 

_JL ' continuing discussions on the nature of the pseudogap, it seems that the preferable "scenario" of its formation is, in 

(— I , our opinion, based on the on the model of strong scattering of current carriers by antiferromagnetic (AFM, SDW) 

Q ' spin fluctuations of short - range order^ [^|3|. In momentum representation this scattering transfers momenta of 

CJ , the order of Q =(-,-) (a - lattice constant of two - dimensional lattice) and leads to the "precursors" of spectrum 

transformation, taking place after the establishment of long - range AFM order (period doubling). As a result we 

T— I , obtain non - Fermi liquid like behavior (dielectrization) of spectral density in the vicinity of the so called "hot spots" 

^ ' on the Fermi surface, appearing at intersections of this surface with "would be" antiferromagnetic Brillouin zone 

^: boundary 0. _ ,_, 

, ■ Within this approach a simplified model of the pseudogap state was actively studied 4, 5], under the assumption 

\^ , that the scattering by dynamic spin fluctuations can be reduced (which is valid for high enough temperatures) to a 
f^ ■ static Gaussian random field (quenched disorder) of pseudogap fluctuations with characteristic scattering vector from 
^T the vicinity of Q, the width of which is determined by the inverse correlation length of short - range order k = £,~'^- 

The review of this approach with applications to the properties of the normal state and for some oversimplified models 
of pseudogap fluctuations influence on superconductivity can be found in Ref. |3] ■ 

In a recent paper 6] we have presented microscopic derivation of Ginzburg - Landau expansion^ and studied the 

influence of pseudogap fluctuations in the model of "hot spots" on the main characteristics of superconducting state 

'"^ (both for s-wave and d-wave pairing), forming "on the background" of these fluctuations. In this paper we have 

^ , considered slightly oversimplified version of our model, where Heisenberg - like spin fluctuations were replaced either 

O ■ by Ising - like, or by spin ~ independent charge (CDW) fluctuations. It was shown that such fluctuations of "dielectric" 

nature generally suppress superconductivity, leading to the drop of superconducting transition temperature and of 

specific heat discontinuity at the transition, as well as to a number of other anomalies. We also discovered two possible 

classes of superconducting order parameter interaction with pseudogap fluctuations, leading to significantly different 

^ ■ scales of suppression of superconducting state. 

5^ i The aim of the present paper is the generalization of the approach of Ref. 0| to the "realistic" case of Heisenberg 

spin fluctuations, as well as the account of (nonmagnetic) impurity (disorder) scattering influence on superconductivity 
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^ Similar charge (CDW) fluctuations also cannot be excluded. 

^ Similar analysis based on Gorkov's equations was given in Ref. Q 



in the pseudogap state. It will be shown that this model allows for semiquantitative modeling of the typical phase 
diagram of a high - temperature superconductor. 

"HOT SPOTS" MODEL AND RECURSION PROCEDURE FOR GREEN'S FUNCTION AND VERTICES. 

Main results of our approach to the "hot spots" model and calculation of a single - particle Green's function in this 
model were described in detail in Refs. [J,l3, methods of calculation for the appropriate vertex parts were discussed in 
Refs. la, Is]- Here we only present the main equations and definitions, giving a short description of changes necessary 
to take into account the spin structure of interaction in Heisenberg model of antiferromagnetic fluctuations. 

Effective interaction of electrons with spin fluctuations in the model of "nearly antiferromagnetic" Fermi - liquid |j| 
is described by dynamic susceptibility, characterized by experimentally determined correlation length ^ and frequency 
ujsf of spin fluctuations, which can depend both on concentration of carriers (and in case of ^ also on temperature), 
Both dynamic susceptibility and effective interaction are peaked (in momentum representation) in the vicinity of 
Q = (tt/o, 7r/a), which leads to the appearance of two types of quasiparticles - "hot" one, with momenta close to the 
points on the Fermi surface, connected by scattering vectors ~ Q, and "cold" one, with the momenta close to the 
parts of the Fermi surface, surrounding diagonals of the Brillouin zone P, Ij, |5{ . 

For high enough temperatures 27rT ^ iVsf and we can neglect spin dynamics |j|. Electron interaction with spin 
(pseudogap) fluctuations reduces then to elastic scattering by appropriate static Gaussian random field. In this model 
it is possible to introduce a simplified form of effective interaction (correlator of this random field) |j, |a|, which 
allows to perform full summation of the Feynman perturbation series, and obtain the following recursion procedure 
determining the single - electron Green's function: 

Gfe(e„p) = ^ . . , ., ^TT 7 (1) 

Sfc(£„p) = W'sik + l)Gfe+i(e„p) (2) 

which is depicted as a symbolic "Dyson like" equation in Fig. ^(a), where we have also introduced: 

iSn - ^k{p) + ikvkK 

Here k = (,~^ is the inverse correlation length of pseudogap fluctuations, £„ = 2TrT{n + 1/2) (for dcfiniteness we 
assume here £„ > 0), 

^ f \vx{p + Q.)\ + \vy{p + Q)\ for odd fc ,, 

^'' ~ 1 |wa:(p)| + \vyip)\ for even k ^ ' 

where v(p) = -4^ is the velocity of a free quasiparticle with the spectrum ^p, which is taken in a standard form |j|: 

^p = —2t{cospxa + cospyu) — At cospxacospyU — /i (6) 

where t is the transfer integral between nearest neighbors, and t' — between second nearest neighbors on the square 
lattice, a is the lattice constant, fi is the chemical potential. 

Parameter W with dimension of energy defines the effective width of the pseudogap. In the model of Heisenberg 
spin fluctuations it can be written as [^]: 

W'^g'^^^g'<in,,^n,,r> (7) 

where g is an interaction constant of electrons with spin fluctuations, < S^ > is the average square of the spin on 
a lattice site, ni-\ , Uii - operators of electron number at a site, with given spin directions. It is clear that similarly 
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(9) 



to correlation length ^, the value of W , in our semiphenomenological approach U,|5|], is also some function of carrier 
concentration (and temperature), to be determined from the experiment. 

The value of s(k) is determined by combinatorics of Feynman diagrams. For the simplest case of commensurate 
charged (CDW) pseudogap fluctuations : 

s{k) = k (8) 

while for the most interesting case of Heisenberg spin (SDW) fluctuations 

^^±^ for odd k 
for even k 

Conditions of applicability of our approximations were discussed in detail in Refs. |j, |5j . 

Remarkable property of our model is the possibility of complete summation of all Feynman diagrams for the vertex 
parts*, describing system response to an arbitrary external perturbation [SJ. Here wejust give appropriate recursion 
relations for "triangular" vertex in Cooper channel, similar to those derived in Ref. ja| in particle ~ particle channel, 
and describing the response to an arbitrary fluctuation of superconducting order parameter (gap): 

A(p,q) = Aqe(p) 

where the symmetry factor, corresponding to the type (symmetry) of the pairing, takes the following form 

1 ( s-wave pairing) 

cospxa — cospyG ( (i^2_j^2-wave pairing) 

and we always assume singlet pairing. The vertex of interest to us can be written as: 

r(e„, -e„, p, -p + q) = rp(£„, -e„, q)e(p) 

Then rp(e„, — e„, q) is determined from recursion procedure of the following form: 



e(p) 



(10) 



(11) 



(12) 



rpfc-i(e„,-£n,q) = 1 ±W^r{k)Gkis„,p + q)Gki-en,p){ 1 



Gj.^(e„,p + q) -G'^\-£„,p) -2ikKVk j 



2ikKVk 



^pk{£n, — £n,q) 



(13) 



which is shown diagrammatically in Fig. n](b). "Physical" vertex is given by rpfc=o(£n, — £n,q)- Additional combi- 
natorial factor r{k) = s{k) for the simplest case of charged (or Ising - like spin) pseudogap fluctuations, analyzed in 
Ref. IqI- For the most interesting case of Heisenberg spin (SDW) fluctuations, analyzed below, this factor is givwn 
by II] (Cf. also Appendix): 



r(fc) 



fc+2 



for even k 
for odd k 



(14) 



The choice of the sign before W^ in the r.h.s of Eq. H13|) depends on the symmetry of superconducting order parameter 
and the type of pseudogap fluctuations [g (see details in the Appendix). Different variants are given in Table I. 

Table I. The choice of the sign in the recursion procedure for the vertex part. 



Pairing 


CDW - fluctuations 


SDW - fluctuations (Ising) 


SDW - fluctuations (Heisenberg) 


s 


+ 


- 


+ 


d 


- 


+ 


- 



It particular, from this Table we can see, that in the most interesting case of d-wave pairing and Heisenberg spin 
fluctuations, we must take this sign "— " , so that alternating signs appear in our recursion procedure for the vertex 
part. However, for the case of s-wave pairing and the same type of pseudogap fluctuations, we must take "-I-", and 
there is no sign alternation in the recursion procedure. It Ref. [g it was shown (for other cases shown in Table I), that 
this difference in the signs behavior of the recursion procedure for the vertex part leads to two cases of qualitatively 
different behavior of all the main characteristics of a superconductor. 



^ Detailed analysis of diagram combinatorics for Heisenberg model of spin fluctuations is given in the Appendix. 
* Including all diagrams with intersecting interaction lines. 



IMPURITY SCATTERING. 

Scattering of electrons by normal (nonmagnetic) impurities can easily be taken into account in self - consistent 
Born approximation, writing down "Dyson" equation for the single - electron Green's function, shown graphically in 
Fig. [3 (a), where in addition to contribution shown in Fig. Q (a), we have included the standard self - energy part 
due to impurity scattering 9]. As a result we obtain the recursion relation for the Green's function as: 

^''^^"^^ " G-i(e„p) - pC/2 ^p G(e„p) - W^sik + l)Gfe+i(e„p) ^^^^ 

where p is concentration of impurities with point - like potential U , and impurity contribution to self - energy contains 
fully dressed Green's function G(e„p) — Gfc=o(£nP)j to be determined self - consistently by our recursion procedure. 
As usual [9|, contribution to this self - energy from real part of the Green's function is reduced to insignificant 
renormalization of the chemical potential, so that l|15|l can be written as: 

^'^^""^^ " ^{en - pU^ Ep /"iG(e„p) + kvun) - Cfc(p) - W^sik + l)Gfc+i(£„p) ^^^^ 

Thus, in comparison with impurity free case, we obtain, in fact, a simple substitution (renormalization): 

£n^ £n- pU"^ ^ /mG(e„p) = EnVe (17) 



1-— E^"^^(^"p) (18) 



If we do not perform fully self - consistent calculation of impurity self - energy, in the simplest approximation we just 
have: 

£n^ £n- PU^ ^ /mGoo(enP) = Snlle = £n + jQSignSn (19) 

p 

where 70 = npU^NQ{0) is the standard Born impurity scattering rate |9( (A'o(O) — density of states of free electrons 
on the Fermi level). 

For "triangular" vertices of interest to us, the recurrence relation with the account of impurity scattering is shown 
graphically in Fig. [3 (b). For the vertex, describing interaction with the fluctuation of superconducting order 
parameter H1Q(1 with d-wave symmetry (jllll this equation is significantly simplified, because the contribution of second 
diagram in the r.h.s. of Fig. El (b) is actually zero, because of validity of J^r, ^(p) = (cf. discussion of similar 
situation in Ref. 10]). Then the recurrence equation for the vertex has the form of Eq. H13|) . where Gfc(±e„p) should 
be taken from Eqs. H15() . (|16() . i.e. just replaced by dressed by impurity scattering, as defined in Fig. 111(a). For the 
vertex describing interaction with fluctuations of the order parameter with s-wave symmetry, we have the following 
equation: 

rpfc_i(e„,-£n,q) = l + /5t/^^G(e„,p + q)G(-e„,p)rp(e, -£„,q) ± (21) 

p 

1 + T^^T, ^ T^^TT ^ W, r Pfc(^"' "^»' °^) 

G^ (£„,p + q)-G^ {-en,P) -2ikKVk j 

where for Gfc(±e„p) we again use Eqs. H15|) . H16|l . and the sign before W^ is determined in accordance with the rules 
formulated above. The difference with the case of the vertex for d-wave symmetry is the appearance of second term 
in the r.h.s. of Eq. (|21|l . which reduces to the substitution: 

l^Vr = l+pU^Yl '^(^"' P + q)G(-£«, P)rp(e, -e„, q) (22) 

p 

Thus our self- consistency procedure can be formulated in this case as follows: start from the "zeroth order" ap- 
proximation G — Goo, Ep = 1, i.e. simply take in Eas. lfT?^ . if^ rje = rjr = I — /of/^/e„ X]p ^"^^oo(£nP)- Perform 



all necessary recursions (starting from some big enough value of k) and determine new values for G = Gk=o and 
Tp = Ffc^o- Calculate again rj^ and r]r using ((T5|l and l|^ . use these values in Eqs. ifTB)) . (|^ etc., until convergence. 
Analyzing the vertex with d-wave symmetry we just take rjr = ^ Sii all stages of calculation. In fact, in this case 
there is no need for achieving full self - consistency over impurity scattering, as it leads only to some small corrections 
to the results of non self ~ consistent calculation using the simplest substitution H19|l 13] • 

SUPERCONDUCTING TRANSITION TEMPERATURE AND GIZBURG - LANDAU COEFFICIENTS. 

Critical temperature of superconducting transition is determined from the equation for Cooper instability of the 
normal phase: 

l-Fx(0;T) = (23) 

where the generalized Cooper susceptibility is defined by diagram shown in Fig. Eland is given by: 

X(q; T) = -TY,Y1 ^(^»P + q)G(-£», ~p)e\p)Tp{s„, -£„, q) (24) 

En P 

Pairing constant V is assumed to be non zero in a layer of the width of 2ujc around the Fermi level, and determines 
the "bare" transition temperature Tcq in the absence of pseudogap fluctuations from the standard BCS equation^: 

l^^^Z^/ '^P- '^PvTT—-2 (25) 

^ „=0 "'0 "'0 ^P ^ " 

where m = [ „ "^ ] is dimensionless cutoff in the sum over Matsubara frequencies. Similar to Ref. 6] all calculations 
were performed for a typical spectrum of bare quasiparticles ^ for different relations between i, t' and /^. E.g. taking 
rather arbitrarily ujc — OAt and Tco = O.Olt it is rather easy to find the value of pairing constant V in (|25|1 . responsible 
for this given value of T^o for different types of pairing. For example, in the case of s-wave pairing we obtain tt- — 1, 
while for the case of dx2_y2--wave pairing we get -^ = 0.55. 

To calculate T^ we need only to know the Cooper susceptibility at q — 0, and this simplifies calculations considerably 
|g|. In general case, e.g. to calculate Ginzburg - Landau coefficients we need x(9; T) at arbitrary (small) q. 

Ginzburg - Landau expansion for the difference of free - energies of superconducting and normal states, written in 
a standard form: 

F, - F„ = A|Aqp + g^CjAqp + ||Aq|4 (26) 

is determined from the usual loop expansion for the free - energy of electrons, moving in a random field of order 
parameter fluctuations (|10|l . 

It is convenient to normalize GL - coefficients A, B, C by their values in the absence of pseudogap fluctuations ja|: 

A^AoKa; C^C'oKc; B = BoKb, (27) 

where 

Ao = No{0)^^ < e2(p) >; Co = ^o(O)^^^ < W{p)fe\p) >; 

Bo = NoiO)^^<e\p)>, (28) 

and angular brackets denote usual averaging over the Fermi surface: < . . . >= j^\q) J2p H^p) ■ ■ •> -^o(O) is again the 
density of state on the Fermi level for free electrons. 



^ We do not discuss the physical nature of this interaction — it may be due to exchange by spin fluctuations, phonons, or by some 
combination of electron phonon and spin fluctuation mechanism. 



Then we have the following general expressions [Q]: 

K, = ^(»^^)-^("^^-) (29) 

Kc = hm x(q;Tc)-x(0;T.) 

i^B = ^ ^ ^ e4(p)(G(£„p)G(-£„, -p))2(rp(£„, -£„, 0))4 (31) 

which were used for direct numerical computations. 

In the presence of impurities all Green's functions and vertices entering these expressions are to be calculated using 
Eqs. lO and (inj. 

The knowledge of GL - coefRcients allows to determine all major characteristics of a superconductor close to the 
transition temperature Tc- E.g. the coherence length is defined as: 

where ^bcs{T) is the value of coherence length in the absence of pseudogap. 
Penetration depth is given by: 

^(^) -( Kb \'' 
\bcs{T) \KaKc) ' ^ ^ 

where we again normalized on the value of \bcs{T) in the absence of pseudogap fluctuations. 
The slope of the upper critical magnetic field close to T^, normalized in a similar way, is given by: 

I <iT It. T, Ka (34^ 



Relative discontinuity of specific heat at the transition is defined as: 

{Cs — Cn)T^o ^cO Kb 

NUMERICAL RESULTS. 

Numerical results for the case of charge (CDW) Ising - like spin (SDW) fluctuations of short - range (pseudogap) 
order were given in Ref. j^. Here we shall concentrate on the analysis of most important and interesting case of 
Heisenberg spin (SDW) fluctuations and the role of impurity scattering (disorder). Due to importance of d-wave 
pairing for copper oxide superconductors, this case is analyzed in more details. 

All calculations (in this section) were performed for typical values of parameters of the "bare" spectrum t' /t = —0.4, 
^/t = —1.3, while for the inverse correlation length we assumed na = 0.2. To spare space we do not present here 
explicit results for dimensionless GL - coefficients Ka, Kb, Kc, giving only the final results for the main physical 
characteristics of a superconductor. 

For dependences on the effective width of the pseudogap we always present these characteristics normalized by the 
appropriate values at T = Tco, while for the dependences on impurity scattering rate 70 we use normalization by 
appropriate values at T = Tco{W), i.e. at the "bare" critical temperature for a given value of W, but in the absence 
of impurities (70 = 0). 



d — wave pairing. 

In Fig. ^we show the dependence of superconducting transition temperature Tc on the effective pseudogap width W 
for several values of impurity scattering rate. It can be seen that pseudogap fluctuations lead to significant suppression 
of superconductivity and for finite disorder we always obtain some "critical" value of W, when Tc becomes zero. This 
suppression of Tc is naturally due to partial "dielectrization" of electronic spectrum in the vicinity of "hot spots" 



Similar dependences for coherence length and penetration depth are shown in Fig. [SJ those for the slope of the upper 
critical field and specific heat discontinuity at the transition — in Fig. El Typically we also observe fast suppression 
of these characteristics by pseudogap fiuctuations. 

Dependences on the value of correlation length of short ~ range order fluctuations is, generally, more smooth — 
in all cases the growth of £, (drop of k) enlarges the effects of pseudogap fluctuations. We do not show appropriate 
results to spare space. 

In Fig. [3 we show dependences of superconducting transition temperature Tc on impurity scattering rate 70 for 
several values of the effective pseudogap width. It is seen that in the presence of pseudogap fiuctuations Tc suppression 
with the growth of disorder is faster than in the absence of the pseudogap {W = 0}, when the dependence of Tc on 70 
in case of d-wave pairing, described by the standard Abrikosov - Gorkov curve [lO|,lll|- Analogous dependences for 
coherence length and penetration depth are shown in Fig. |S1 while for the slope of the upper critical field Hc2{T) and 
specific heat discontinuity — in Fig. |51 Again we see, that impurity scattering (disorder) leads to fast suppression of 
the last two characteristics, e.g. the effect of pseudogap fluctuations is enhanced by impurity scattering. 

Dependences on parameters of the pseudogap, obtained above, are qualitatively similar to those found in Ref. [a| 
for the case of charged (CDW) pseudogap fluctuations, when we also are dealing with recursion procedure for the 
vertex part with alternating signs. At the same time certain differences appear here due another combinatorics of 
diagrams. Dependences on the impurity scattering rate (disorder) in this model were not studied previously^. 

These dependences are in qualitative agreement with majority of the data, obtained in experiments on copper 
oxide superconductors in the pseudogap region of the phase diagram (from underdoped to optimally doped region) . 
Below we shall demonstrate that these results may be used for direct modeling of the typical phase diagram of these 
compounds. 

s — wave pairing. 

Analysis of the case of s-wave pairing is interesting mainly as demonstration of characteristic differences with the 
case of d-wave pairing. At present, there are no definite experimental data on s-wave superconductivity in systems 
with pseudogap, though it is quite possible that such systems will be discovered in some future. 

Our calculations show that pseudogap fiuctuations significantly suppress superconductivity also in this case (Fig. 
I10|l . though the characteristic scale of these fluctuations, necessary for a major suppression of the superconducting 
state, is much larger here than in the case of d-wave pairing. This result was obtained earlier in Ref. pj, though 
it is necessary to note that in the case of Heisenberg (SDW) fluctuations, considered here, there is no characteristic 
"plateau" in the dependence of Tc on VF, which was obtained for the case of charged (CDW) pseudogap fiuctuations 
in Ref. [y|- On the same scale of W takes place also the major suppression of specific heat discontinuity at supercon- 
ducting transition, as shown at the insert in Fig. EH Similar dependences for the coherence length and penetration 
depth are analogous to those obtained in Ref. 6] and we do not show these here. In FigE]we show dependences of 
Tc on impurity scattering (disorder). Besides relatively weak suppression of Tc by disorder, connected [7| with small 
suppression of the density of states at the Fermi level by disorder, we can observe also a weak growth of Tc with 70, 
which is apparently due some "smearing" of the pseudogap in the density of states. 

In Fig. El we show impurity effects for the case of s-wave pairing on the coherence length and penetration depth. 

Finally, in Fig. El we demonstrate the influence of impurity scattering (disorder) on the slope of the upper critical 
fleld and specific heat discontinuity. Again the specific heat discontinuity is significantly suppressed by disorder, while 
the slope of Hc2{T) demonstrates qualitatively different behavior, than in the case of the d-wave pairing: the growth of 
disorder leads to the growth of this slope, as in the standard theory of "dirty" superconductors [23, while pseudogap 



Analysis of impurity dependence of T^ in Ref. JTl was performed for the case of recursion procedure without sign alternation, appearing 
for the case of Ising — like SDW fiuctuations, when superconductivity suppression is much slower. 
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fluctuations further increase the slope of Hc2(T). In the absence of pseudogap fluctuations, similar differences (between 
s-wave and d-wave cases) in the behavior of the slope of Hc2{T) under disordering were noted earlier in Ref. |lO|. 

MODELING OF THE PHASE DIAGRAM. 

The described above model of the influence of pseudogap fluctuations on superconductivity allows to perform 
simple modeling of a typical phase diagram of high - temperature superconducting cuprates^. First attempt of 
such modeling in an oversimplified version of our model was made in Ref. |l2l |. The main idea is in identifying 
our parameter W with an experimentally observable effective width of the pseudogap (crossover temperature into 
pseudogap region of the phase diagram) Eg ^ T* , which is determined by numerous experiments [ll Ij, y{ - It is well 
known that this crossover temperature is practically linear in concentration of doping impurity (carrier concentration) , 
starting at small concentration from the values of the order of lO^K and going to zero at some critical concentration 
Xc ~ 0.19. .0.22, which is slightly higher than the "optimal" value Xo ~ 0.15. .0.17 [llll2|- Accordingly, we can assume 
similar concentration dependence of our pseudogap parameter W{x)^. In this sense we can say that concentration 
dependence W{x) directly follows from the experiment. Then the only parameter of the model to be determined 
remains the concentration dependence of the "bare" superconducting transition temperature Tcq{x), which would 
have existed in the absence of pseudogap fluctuations. The knowledge of this dependence allows then to calculate the 
observed concentration dependence of superconducting transition temperature Tc{x) solving equations of our model. 
Unfortunately, as was already noted in Ref. |a|, concentration dependence Tco{x) is, in general case, unknown and is 
not determined from any known experiment. Thus, it remains the fitting parameter of our theory. 

Assuming following Ref. 12] that Tco{x) can also be described by a linear function of concentration x, going to 
zero at a; = 0.3, and varying the value of Tco{x = 0) to obtain the correct value of Tc{x = Xo) it is possible to calculate 
the "observed" dependence of Tc{x). An example of such calculation for the case of d-wave pairing and scattering 
by charged (CDW) pseudogap fluctuations [a|, using typical dependence of W{x), is shown in Fig. ^J It can be 
seen that even with such crude and rather arbitrary assumptions, the "hot spots" model allows to obtain the Tc{x) 
dependence, which is pretty close to the experimentally observed. Similar calculations for the case of Ising - like 
model of interaction of electrons with spin fluctuations (no sign alternation in the recursion procedure for the vertex 
|a|) show that the reasonable values of Tc{x) can only be obtained for rather unrealistic values of W{x) about an 
order of magnitude higher than those observed experimentally. 

Within our BCS - like model for the "bare" Tco, an assumption of significant concentration dependence of this 
parameter seems rather unrealistic^. Thus, let us assume that the value of Tco has no dependence on carrier (doping 
impurity) concentration x at all, but take into account the fact that introduction of doping impurity immediately 
leads to the appearance of some random impurity scattering of carriers (due to internal disorder) , which can be simply 
described by introduction of some linearly growing impurity scattering rate 7(x). Now let us just assume that this 
growth of internal disorder leads to complete suppression of d-wave pairing at a; = 0.3, in accordance with Abrikosov - 
Gorkov dependence [lOllllI- Results of calculations for the phase diagram oi La2~xSrxCuOi ~ type system within this 
model, and for the case of Heisenberg pseudogap fluctuations and the described role of impurity scattering is shown 
in Fig. ^] The values of different parameters used in this calculation are shown at the same figure. "Experimental" 
values of Tc (a;), shown in this figure (as well as in Fig. I14|l by "diamonds", were obtained from the empirical relation 
[HEi: 

^'^""^ T - 1 - 82.6(a; - x^f (36) 



Tc{x = 

which represents rather good fit to experimental data for concentration dependence of Tc in a number of cuprates. We 
can see that in all underdoped region our model gives practically ideal description of "experimental" data for quite 
reasonable values of W{x). At the end of overdoped region agreement is less satisfactory, but it should be noted that 
both the relation (|36(l is not quite good here, and our model of superconductivity suppression in overdoped region 



"^ We "neglect" here the existence of a narrow region of antiferromagnetic ordering in Mott dielectric appearing at low concentrations of 

doping impurity, limiting our analysis to a wide region of a "bad" metal. 
* Obviously, this identification can be done up to some unknown factor ~ 1. 
^ Within this approach any dependence of TcO on x can only be due to some weak concentration dependence of the density of states on 

the Fermi level. 



is obviously too crude. Also it should be noted that we have not performed any special fitting of parameters of our 
model to increase agreement in this region. 

It is interesting to analyze the behavior of superconducting transition temperature T^ under additional disordering at 
different compositions (concentrations of carriers). There were numerous experimental studies with such disordering 
being introduced by additional impurities [la, ll3l or by irradiation by fast neutrons [13 and electrons [13, Il9l| . 
Special discussion of the role of such additional disordering in the context of the existence of the pseudogap state was 
presented, apparently, only in Ref. |l7l |. 

In our model such additional disordering is easily simulated by the introduction of an additional impurity scattering 
rate 70, which is just added to our parameter of an internal disorder ^{x). Results of our calculations for two values of 
this parameter are also shown in Fig. 1151 We can see, that in complete agreement with experiments |l7j . introduction 
of additional "impurities" (disorder) leads to rather fast reduction of the superconductivity region on the phase 
diagram Also in complete agreement with our conclusion made above with respect to Fig. |7| and with the results of 
experiments |1!tI Il7l | , superconductivity suppression by disordering in underdoped (pseudogap) region is significantly 
faster, than for optimal doping. It could be expected that introduction of "normal" disorder, obviously leading to 
some suppression of the pseudogap in the density of states, could have lead to some "slowing down" of Tc suppression. 
However, our calculations show that such effect in the case of rf-wave pairing is just absent. 

The problem, however, is that in all cases our calculations show that Tr suppression is faster than it is predicted 
by the standard Abrikosov - Gorkov curve in case of d-wave pairing |lClllll |. At the same time, most attempts to fit 
experimental data on disordering in cuprates |lfiL llSl Il9l | , show that such suppression is in fact significantly slower, 
than predicted by Abrikosov - Gorkov dependence. This fact remains one of the major unsolved problems of the 
theory of high - temperature superconductors |2i| . One of the possible solutions of this problem may be connected 
with rigorous description of disordering effects in superconductors, which belong to transition region from "large" pairs 
of BCS theory to "compact" pairs (bosons), appearing in the limit of very strong coupling |2l|. Another interesting 
possibility to explain the "slowing down" of Tc suppression is connected with an account of anisotropy of impurity 
scattering, analyzed in detail in Refs. [l(j,|22|. This effect can be rather easily included in our calculational scheme. 
It seems especially significant in connection with the established strong anisotropy of elastic scattering (with d-wave 
symmetry), observed in ARPES experiments on Bi2Sr2CaC'u20s+s |23,|2J|. Appropriate scattering rate changes in 
the interval of 20 — 607716^ 24], which is nearly an order of magnitude higher, than the maximal value of 7(2:), used 
in our calculations, which once again demonstrates an unusual stability of d-wave pairing in cuprates towards static 
disorder. It should be noted that our model of electron self - energy in fact describes analogous anisotropy of elastic 
scattering, corresponding to its growth in the vicinity of "hot spots" , but we do not observe the effect of the "slowing 
down" of Tc suppression in our calculations. 

Our results show that despite the obvious crudeness of our assumptions, the "hot spots" model allows to obtain 
rather reasonable (even semiquatitative) description of the region of existence of superconducting state on the phase 
diagram of high - temperature superconducting cuprates^*^. The main deficiency of our approach remains rather 
indeterminate scenario for concentration dependence of the "bare" superconducting transition temperature. 

CONCLUSION 

The analysis presented above shows, that the model of the pseudogap state, based on the concept of "hot spots" , 
can provide a reasonable description of the main properties of the superconducting phase of cuprates, as well as 
of the phase diagram, using relatively small number of fitting parameters, most of which can be determined from 
independent experiments. 

Let us stress that all calculations were performed using the standard assumption of self - averaging nature of 
superconducting order parameter (gap) over the random fields of impurities and pseudogap fluctuations. Usually this 
assumption is justified for superconductors with coherence lengths (size of Cooper pairs) much exceeding the other 
microscopic lengths in the system, such as mean free path or correlation length of pseudogap fluctuations ^. 

In the class of models of the pseudogap state under consideration this assumption is not, in general, valid and 
important effects due to non self - averaging |i25t i2fi| | , which lead to a qualitatively new picture of inhomogeneous 
superconducting state, with "drops" of superconducting phase existing for temperatures T > Tc- In principle, there are 



^^ A bove we always assumed that we are dealing with hole - doped cuprates, where concentration dependence of T* (x) is well established 
Q II3I . For electron doped systems the data on the pseudogap state are rather fragmentary. 
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experimental data, directly supporting such picture of inhomogeneous superconductivity in cuprates [23, 1281 |29j . Of 
course, we are far from claiming, that these experiments directly support oversimplified theoretical models developed 
in Refs. [25|,|26j. At the same time, these results stress the importance of rigorous analysis of non self - averaging 
effects in relatively "realistic" models of the pseudogap state, such as the "hot spots" model, described above^^. 

This work was supported in part by the RFBR grant No. 02-02-16031, as well as by the program of fundamental 
research "Quantum macrophysics" of the Presidium of RAS, the program "Strongly correlated electrons in semicon- 
ductors, metals, superconductors and magnetic materials" of the RAS Department of Physical Sciences, and also 
under the project of the Ministry of Education and Science of the Russian Federation "Studies of collective and 
quantum effects in condensed matter" . 



^^ The picture of "superconducting drops", existing above Tc allows to understand a number of experiments, which are usually interpreted 
as supporting the superconducting nature of the pseudogap in cuprates. 
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DIAGRAM COMBINATORICS FOR HEISENBERG PSEUDOGAP FLUCTUATIONS. 

To analyze diagram combinatorics, let us consider the limit of infinite correlation length of pseudogap fluctuations. 
In this case the spin density scattering electrons takes the form: 

Sq = S^(q-Q) (37) 

and averaging over Gaussian spin fluctuations reduces to the usual integration |j] : 

<...>= ^, / dSe-^^ . . . (38) 

Accordingly, in this limit we can first solve the problem for an electron, moving in the coherent field of spin density 
(|37|l . and afterwards perform averaging (|38|l over its fluctuations. For further analysis it is convenient to introduce 
fluctuating field 5 — -y=S — "potential", scattering electrons. Then averaging (|38() over spin fluctuations reduces to 
the averaging over fluctuations of this random field: 



3 ^ 


J. 


/ dSie" 




27rM/2^ 


2nW^J, '^^io ' " ' 



We see that actually we have two fluctuating fields scattering electrons: the real longitudinal field 5i = -4=8^ and 
complex transverse St, characterized by the amplitude \St\ and phase ip, connected with two transverse components 
of S. 

Such averaging generates diagram technique with two types of effective interactions y. One can be represented 
say by dashed interaction line and correspond to: 

Veffi = y < ^.q^.-q >- ±— ^(q - Q) (40) 

where the sign "— " is taken when spin projection under this line is changed (e.g. when this dashed line surrounds an 
odd number of spin - flipping operators 5+, S-), and another, represented by wavy line: 

Veff2 = y < 5+q5— q >= 2^<5(q - Q) (41) 

The averages like < 5*+ 5*+ > and < S-S- > are zero after averaging over the phase in H39|l . 

Let us start with electron moving in coherent flled of spin density H37() . In this case single - particle Green's function 
is represented by 2 x 2 matrix with 4 independent components"'^^, which are determined by the following system: 

Gii-11 = Gi + Gi(5;G2-f;i-f + GiStG2i-ii 

G2|;i| — G2(5;Gif ;i| + G2(5tGij^;il 

G2j.;it — — G2(5/Gi4;i| + G2i5^Gi|;i| 

Gi|;i| = — Gi^/G2X;i| + Gi(5j*G2t;i| (42) 

where we have introduced short notations: (e„,p) -^ 1, {£n,P + Q) ^ 2 and Gi = jt-^it-, G2 ~ ^^ \ . Then we 

obtain: 

Gl|:l| 

^l'^2'-\0r Uj^-U^- -\0\' 

A* 

(43) 



G2' 


r Si 




Gr^G2-i-|<5P' 
Gi;;it — 0; 


^^'^^' Gr^G2-^- 


-\s? 


^^^=^' G^'G^'- 


~\S\^ 



where \5\ = y^Sf + |<5tp is the amplitude of the fleld 6. 



^^ Components diflFering from these by the change of all spin projections can be obtained by substitution: <5; — ► —(5; , St 
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The averaged Green's function is given now by: 

-ITITT^ / d\5\\5?e-'-^ ^,-,^\ „„ (44) 

This integral form can be easily represented [J[ by continuous fraction Q, |5J with k = with combinatorial 
coefficients s(fc), defined in ^. 

Situation with combinatorial factors r(fc) for two - particle vertices is more complicated. In fact we have to consider 
four types of vertices, shown in Fig. 1161 For all these vertices the recursion procedure has the form of 1)13(1 . but signs 
and combinatorial factors r(fc) are different. Consider first all these vertices for an electron in coherent field 5. 

1. Charged vertex (spin projection is conserved) in diffusion (particle - hole) channel (Fig. llGlfal): 

j.ch_\^r' r - ('^2^2')'^ + W i.r\ 

id — 2^Lrl|;jo-Lri'cr;l'T — ^ (45] 

where i and a take the values 1, 2 and t, ii and we introduced notation (£'j,p') -^ 1', (eJi,p' + Q) ^ 2' ds = 

2. Charged vertex in Cooper (particle - particle) channel^'^ (Fig. 1161 fb')'): 

Tc = / ^Gi|:icrGi'|:i'o- = (46) 

as 

3. Spin vertex (spin projection changes sign) in diffusion (particle - hole) channel (Fig. 1161 ic)'): 

r'^p -\^ r r - (^2^2')" - 5i 

i- d — 2^tjlT;io-'-^i'-(T;l'X — -j [^1) 

Z,fT 

4. Spin vertex in Cooper (particle - particle) channel (Fig. 1161 fd)): 

rsp sr^ r< n (^^2^2')'^ + (|^t|^ - ^; ) , . 

^ c — 2_^'^\\\i'y'^Vi\i'-<y — -, i48j 

2,0" 

Physical vertices are obtained from these by averaging by fluctuations of coherent field (5 with the help of ((39(1 . 
Finally, we can see that the vertex Y^J^ is determined by Eq. ((45(1 . while all other vertices take the form^^: 

,.JE22iri^ ,49) 

where "+" corresponds to vertices TJ;'' and F^^, while "— " is taken for F^''. 

In case of the vertex FJ^'' we have, obviously, r(fc) — s{k). It follows directly from the fact that diagram expansion 
for the physical vertex < F^^ > can be obtained by inserting of appropriate bare vertex into all electron lines in any 
arbitrary diagram for the single - particle Green's function. This insertion docs not change the direction of either 
electronic line or spin projection, so diagram combinatorics does not change. 

In the limit of infinite correlation length any "skeleton" diagram for the vertex differs from "ladder" diagram of the 
same order in interaction ■^(5(q — Q) only by sign and the factor of 2^, where p — is the number of wavy interaction 
lines. Thus, the sum of all "skeleton" diagrams in the given order reduces to the contribution of the appropriate 
"ladder" diagram with interaction ^^(5(q — Q), multiplied by combinatorial factor, which we call the "number" of 
"skeleton" diagrams of the given order. 

First term in Eqs. ((45() - ((48(l is the same for all vertices and generates, after averaging procedure, the "numbers" 
of "skeleton" diagrams of even order in W^ (this term corresponds to i = 1 contribution in these expressions) . We 
see that the "numbers" of "skeleton" diagrams of even order are the same for all types of vertices. The second 



^^ It appears in case triplet pairing. 

^"^ This form is equivalent to I4bi . 1471 . 1481 after averaging 
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term in these expressions generates "numbers" of diagrams of odd order (corresponding to contributions with i = 2). 
Accordingly, the "numbers" of "skeleton" diagrams of odd order for three vertices, defined by (|^ . is equal to ±i of 
appropriate "numbers" for the vertex F^''. The sign "— ", corresponding to the vertex F^^, may be compensated by 
appropriate change of sign in recursion procedure for this vertex. It follows that the sign before the second term in 
(|49|l determines the sign in recursion procedure H13|) for these vertices, while combinatorial factors r(fc) for all three 
vertices are the same. 

The "number" of "skeleton" diagrams of the order L is now^^: 

3^ n '^^) (50) 



l<k<L 



Thus, for even L — 2n we obtain: 



For odd L ^ 2n + 1: 



n rik)= H sik) (51) 



l<k<2n l<k<2n 



n ^w=^ n ^(^) (52) 



3 

l<k<2n+l l<k<2n 

Then, taking into account 0, we immediately obtain H14(l . 

In this paper we were mainly interested in the vertex F^^. The analysis, performed above, shows that for this vertex 
we get the recursion procedure without alternating signs for the case of s-wave pairing, when symmetry factor e(p) 
of the vertex is equal to 1. In case of d-wave pairing, when superconducting gap changes sign after change of the 
momentum by Q (i.e. e(p) = — e(p + Q))), the sign in recursion procedure is changed to the opposite |6j and we 
obtain recursion procedure with alternating signs. Note, that in the case of Ising - like spin fluctuations, analyzed in 
Ref. |a|, situation with signs in recursion procedure for the vertex is strictly opposite. This fact is easily understood 
from an expression H48|l for the vertex F^^. In Ising case two transverse components (i.e. the field St) just disappear, 
leading to the sign change of the second term in H48|l and in recursion procedure. 



^^ The factor of 3^ appears due to the fact that both recursion procedure 1131 and combinatorial coefficients r{k) correspond to expansion 
in powers of W^, while the "number" of "skeleton" diagrams is defined for an expansion in powers of ^^- 
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FIG. 1: Recurrence equations for Green's function (a) and "triangular" vertex (b). 
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FIG. 2; Recurrence equations for Green's function (a) and "triangular" vertex (b) including scattering by random impurities. 
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FIG. 3: Diagrammatic representation of the generalized Cooper channel susceptibility x{^) 
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FIG. 4: Tc. dependence on the effective width of the pseudogap W for the case of d-wave pairing and for several values of 
impurity scattering rate 'jo/Tco: 0-1; 0.18 - 2; 0.64 - 3. Inverse correlation length of short - range order Ka=0.2 
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FIG. 5: Dependence of the square of coherence length and penetration depth on the effective width of the pseudogap W for 
the case of d-wave pairing and for several values of impurity scattering rate yo/Tco- 0-1; 0.18 - 2; 0.64 - 3. 
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FIG. 6: Dependence of the slope of Hc2{T) and specific heat discontinuity at superconducting transition on the effective width 
of the pseudogap W for the case of d-wave pairing and for several values of impurity scattering rate ^o/Tco: 0-1; 0.18 - 2; 
0.64 - 3. 
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FIG. 7: Tc dependence on the impurity scattering rate (disorder) 70 for the case of d-wave pairing and for several values of the 
effective pseudogap width W/Tco: 0-1; 2.8 - 2; 5.5 - 3. 
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FIG. 8: Dependence of the square of the coherence length and penetration depth on the impurity scattering rate (disorder) 70 
for the case of d-wave pairing and for several values of the effective pseudogap width W/Tco: 0-1; 2.8 - 2; 5.5 - 3. 
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FIG. 9: Dependence of the slope of Hc2{T) and specific heat discontinuity at superconducting transition on the impurity 
scattering rate (disorder) 70 for the case of d-wave pairing and for several values of the effective pseudogap width W/Tco- - 
1; 2.8 -2; 5.5-3. 
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FIG. 10: Tc dependence on the eflective width of the pseudogap W for the case of s-wave pairing and for two values of the 
impurity scattering rate ^o/Tcq: 0-1; 20-2. Inverse correlation length of short - range order Ka=0.2. At the insert — 
characteristic dependence of specific heat discontinuity for the same parameters. 
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FIG. 11: Tc dependence on the impurity scattering rate (disorder) 70 for the case of s-wave pairing and for several values of 
effective pseudogap width W/Tco: - 1; 8 - 2; 15 - 3. 
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FIG. 12: Dependence of the square of coherence length and penetration depth on the impurity scattering rate (disorder) 70 for 
the case of s-wave pairing for two values of the effective pseudogap width W/Tco- -- 1; 15 - 2. 
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FIG. 13: Dependence of the slope of Hc2{T) and specific heat discontinuity at the transition on the impurity scattering rate 
(disorder) 70 for the case of s-wave pairing and for two values of the efi'ective pseudogap width W/Tco: 0-1; 15-2. 
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FIG. 14: Model phase diagram for the case of pseudogap fluctuatfons of CDW type (d-wave pairing) and "bare" superconducting 
transition temperature Tco with linear dependence on carrier concentration. 
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FIG. 15: Model phase diagram for the case of Heisenberg (SDW) pseudogap fluctuations (rf-wave pairing) and "bare" super- 
conducting transition temperature Tco independent of carrier concentration, taking into account the role of internal disorder, 
linear over the concentration of doping impurity 7(x)- 
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FIG. 16: Vertex parts with different diagram combinatorics. 



